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We examine spin-Peierls systems in the presence of spin defects which are introduced by 
replacing magnetic ions Cu^^ with non-magnetic ones Zn^"*" in CuGeOs. By using the action 
for the bosonized Hamiltonian, it is shown directly that the antiferromagnetic state induced 
by the spin defects coexists with the spin-Peierls states. Further the doping dependences of 
both transition temperature of spin-Pcierls state and the spin gap have been calculated. The 
transition temperature of the present estimation shows good agreement quantitatively with 
that observed in Cui-^Zn^Oa for the region of the doping rate, 5 < 0.02. 
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§1. Introduction 

The discovery of inorganic spin-Peierls (SP) system CuGeOs, which has chains of Cu^^ with 
spin 1/2 along c-axiJ3), has attracted much attention to the role of disorder in SP systems. In the 
material, the antiferromagnetic (AF) exchange energy between nearest neighbor spins is estimated 
as J = 183K from the magnetic field which corresponds to saturation of magnetization!^. The 
compound in the absence of disorder undergoes SP transition at Tsp = UlS , below which the 
lattices spontaneously dimerize and the magnetic susceptibility rapidly drops to zero along all the 
axes owing to appearance of the energy gap between the singlet ground state and the triplet excited 
one. The spin gap is estimated as = 2.1meV from the inelastic neutron scattering^. 

The effects of the disorder on SP systems have been mainly investigated on the materials where 
Cu ions are replaced by ZnBHilBillHlllst or Ge ions are replaced by Si0'0'liil'. It has been observed 
that the doping decreases Tsp. Further, the fact that a new AF ordered phase appears below Tsp 
has been now well established from measurements of the magnetic susceptibility and the neutron 
scattering. This fact shows that the AF states coexists with the SP state. 

* E-mail: h44770a@nucc.cc.nagoya-u.ac.jp 
** E-mail: e43428a@nucc.cc.nagoya-u.ac.jp 
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In terms of the phase representation of SP Hamiltonian, Fukuyama et al. exammed the states 
between the two impurities by assuming the boundary condition that the lattice dimerization is 
suppressed at the location of the impurities0\ Treating classical equations with renormalized 
quantum fluctuation, they have shown that the AF state coexists with the SP state in the presence 
of disorder. 

In the present paper, we investigate the effect of the spin defects on SP states without the 
assumption for the boundary condition by deriving the action for the SP system and show the 
coexistence of the two kinds of states. It turns out that our results are consistent with those given 
by Fukuyama et al. . Further, we examine the variation of both Tsp and A5 by the increase of 
the doping rate. In §2, we derive the action of the SP systems with the spin defects from Peierls- 
Hubbard model. The transition temperature is calculated in §3. In §4, we prove the coexistence 
of AF state and SP one and investigate the doping dependence of A5. Section 5 is devoted to 
summary and discussion. 

§2. Model 

Since the replacement of Cu with Zn gives rise to the variation of both the charge and spin degrees 
of freedom, we examine the effects of the spin defects by starting the Peierls-Hubbard model in 
one-dimension at half filling instead of the conventional SP model. The Hamiltonian is given as 



n = -t'^{l + X{ui - Ui+i)} + h.c. 



+ ^Y1 4,TCi,T4,iC«a + ^ ^i^l - ui+if, (1) 

where t and U are the hopping and the on-site Coulomb energy, respectively. The quantity A is 
the coupling constant of the electron-lattice interaction and c]^^ is a creation operator of electrons 
with the site / and the spin a. The quantity K is the elastic constant with ui being the lattice 
distortion at the site I. 
The model in eq.(l) can be expressed by the phase variables as followsEl\ 

Vn I . f 1_ 



+ / dxcos20^ 



+ ^'^^'^ I dxcos; 



dxu(x) sin^+ cos</>+ H / dxn^(x), (2) 

where a is the lattice constant, B cx X and u{x) = {—lYui with x = la. The quantity a denotes 
the cut-off for the large wavenumber. The velocity of excitation of the charge (spin) is given by 
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(^cr)- The quantities Qp and g^^ are the coupling constants of the interaction for the umklapp 
scattering and the backward scattering respectively. The coefficients Kp and K^j which also include 
parameters of interaction characterize one-dimensional system. The phase variables, 9± and <j3±, 
express fluctuations of the charge and the spin degrees of freedom, respectively, and satisfy the 
commutation relations, [9^{x),0^{x)] = [(f)^{x),(j)-{x)] = i7rsgn(x — x') and the others are zero. 
Note that we assumed t > U in deriving eq.(2). The density of the z component of the spin is 
expressed by dx4>+/ {^tt) — (— 1)' sin0+ sin(/)_(_/7ra, where the first and the second term express the 
slowly varying part, S^, and the staggered one, 5J, respectively. 

In the above phase Hamiltonian, the charge excitation has a gap and the phase variable 0+ is 
fixed to the value of 0+ = (n-|- l/2)7r ( n : integer ). Thus, eq.(2) in the absence of disorder becomes 
the same as that of phase representation of the conventional SP Hamiltonian©. The fact indicates 
that the regime of the weak interaction is analytically connected to that of the strong interaction. 
Therefore the parameters in eq.(2) can be read as those of phase Hamiltonian of SP model. 

The defects of spins are introduced as follows. At the location where Cu^'^ are replaced by Zn^^, 
there appears the charge defect which is expressed by as a kink of 0+. Since the fluctuation of the 
charge density is given by dxOj^/-K, the quantity 0+ jumps by vr 0) at the defect where sin 0+ varies 
from ±1 to =f1 resulting in the change of the sign of the coupling between spin and lattice. Thus 
the action of SP systems in the presence of the spin defects is obtained as 



|- j drdx|(a,0+)2 + -L(a.0^ 



AttK, 

H — / drdx cos 2(j)j^ 
vra^ J 

drdxm(2;)ti(x) cos 0+ + /3 / dxu^(x), (3) 

where /? is the inverse of temperature, T. The quantity m{x), which satisfles |m(x)| = 1, 
changes the sign at the location of the spin defects. The staggered component of S'^ is given 
by — (— l)^m(3;) sin(/)+(x)/7ra. In the following, we study the case where the spin systems are ex- 
pressed by Heisenberg chains, where fg- = 7rJa/2lll^, K^j 1 due to the rotational symmetry in 
spin spaceEl^ and 

§3. Temperature of Spin-Peierls Transition 

We calculate Tsp from the softening of phononS^ with q = tt. The complete softening takes 
place at T = Tsp resulting in the vanishing of the effective elastic constant, K^jf, which comes 
from the renormalization of K in eq.(3) through the electron-phonon coupling. By integrated out 
the spin degree of freedom up to the second order, K^ff is given as follows, 

ft 

Keff = K- —B^ / da;/i(x)e-^/«i^(e-2^'/«), (4) 

UJc Jo 
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where ujc = Va-a"^, = v^/{itT) and K{k) is the complete ehiptic integral of the first kind. The 
function h{x) is defined by (l/L) J dX {m{X + x/2)m{X — with (• • expressing the 

impurity average. By noting that h{x) is given by l — 26\x\/ a for 25\x\/a <^ 1 and h{x — > ±00) = 0, 
we make use of the approximation given by e"^'''^'/" where 6 denotes the doping rate. Then 
Tsp/Tgp with Tgp being the transition temperature in the absence of the impurities is determined 
self-consistently by the following equation, 

^ = 1 duu^Jmp)mp/np)K{u\ (5) 
Tgp Ci Jo 

where Ci = Q dtK{t'^) ~ 1.72. The asymptotic behavior is given as follows, 

Tsp 

for 

2C,^--,]\ tor ^»1 • 




-'SP /J ''-SP 



(6) 



where C2 = — Jq dt In tK{t'^) ~ 1.59 and 7 ~ 0.577 is Euler's constant. We show Tsp/T^p as 
a function of J5/Tgp in Fig. 1(a). According to the measurements of magnetic susceptibility 
and neutron scattering, Tsp of Cui-^Zn^GeOa decreases linearly for 5 < 0.02 and saturates from 
6 ~ 0.02 to at least 6 ~ O.O5i'0'. In Fig.l(b), by substituting the values of J = 183K and 
T^p = lAK into eqs.(^) and (6), the present result (solid curve and dashed line) is compared with 
the experimental one given by the measurements of the neutron scattering (open circles) and the 



magnetic susceptibility (filled circles) which is shown in Fig.l of ref.lC. Figure 1(b) shows good 
agreement between the present result and the experimental one for (5<0.02. However, our result 
cannot explain the saturation of the transition temperature in the region of the doping rate of 
0.02 <6 < 0.05. 

§4. The Case ofT <^Tsp 

In this section, we examine the coexistence of the AF state and the SP state in the presence of 
the spin defects and calculate the doping dependence of A5 based on the action given by eq.(3). 

We note that the change of the sign of the spin-lattice interaction at the location of the 
spin defects can be reexpressed as m(x)cos04.(x) cos{(p+{x) + n(x)7r} and m{x)sm(p+{x) — > 
sin {i;^>+(x) + n(x)7r} where the integer n{x) jumps alternatively by ±1 at the spin defects. By 
rewriting (f)^ + n{x)TT as (p, the action, S, can be given as follows, 

\2 
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Fig. 1. (a) The critical temperature of spin-Peierls transition normalized by that without doping as a function of 
{J/T^p)5. The dotted (dashed) curve shows an asymptotic form for JS/Tgp <g 1 {';^ 1) given by eq.(6). (b) 
Comparison of the present result with the experimental result (open and filled circles) in ref.^ where J — ISS-K" 
and Tgp — 14K. Here the open and the filled circles expresses Tsp determined by the measurements of the neutron 
scattering and the magnetic susceptibility, respectively. 



where Xi is the position of the spin defects. With the use of the number pi = ±1, dxn{x) = 
J2iPi^ix ~ X'l) and PiPi+i = —1. In eq.(7), we divide the lattice distortion as u{x) = —u + u{x) 




(7) 
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where u{x) denotes the spatial variation around the uniform value, u. Correspondingly, the action 
S is divided as S" = Si[u] + S2[u,u{x)]. We apply the variational method to 5'i[n], whose trial 
action So[u] is given by 



So[u] = ^J drdx + ^oV' + ^J9r<Pf] 



Here L is the system size, and qo is related to the spin gap in the presence of the disorder A5 as 
As = Vfjqo- The quantities, go and u are determined by minimizing the following Free energy©, 

Ftr = Fo-Bujdx (cos 4>)so - ^ / ' 

where (■ ■ ■) expresses the average with respect to 5*0 [u]. On the other hand, u{x) is determined 
by minimizing the quantity, — In (exp { — S'2[n, n(a;)]})g^ with respect to u{x). 

In order to calculate the various quantities averaged by ^o, we first calculate the generating func- 
tion, l^[J(x,r)] = (exp{— i / drdx J(2;, r)(;/)(x, r)})^^ with J(x,r) being an arbitrary real function. 
By applying path integral method, the generating function is given as follows, 

W[J{x,t)] = exp I ~ drdxdr'dx' 

H{x -x',T- t'; T)J{x, r) J(x', r')| 
X exp|-i J drdxF(x)J(x,r)| , (10) 



where 



H{x — x',T — r'; T) 
1 exp{-ig(x - xQ + ia;n(r - tQ} 

q,^r. —{q +qQ+u;n ) 



and 



^ vr igexp {ig(xj — x)} 



= |E^'ie-'^°'"""''sgn(x-xO, (12) 

i 

with ujn = 2mrT ( n : integer ). Equations (10)-(12) lead the following results, 

S:{x) = -t^ {sin <t>)s. 
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'^^'exp|-iii'(0,0;r)|sinF(x), (13) 

= -|f exp|-^i7(0,0;r)|cosF(x). (14) 

with X = la. Here the first hne in eq.(14) is derived by minimizing — In (exp {— S'2[n, u(x)]})^^ 
with respect to u{x). The function |-F(x)| takes almost a value of tt/2 near the location of the 
spin defects and decreases rapidly to zero far from the location. Thus, without the assumption of 
the boundary condition, we obtain that the lattice dimerization is suppressed and AF moments 
appear near the location of the spin defects. The characteristic length of the spatial dependence 
is given by 1/qo- Equations (13) and (14) show the facts that the phase (j){x) can be divided 
as (j){x) = 4>ci{x) + 4>q{x) where (pdi^) is the classical part and (pq{x) expresses the quantum 
fluctuation around 4>ci{x) where (pdix) = F{x) and (^(j)g{x)'^ = H{0, 0; T). In Fig.2, we compare the 
classical part of the phase obtained in the present study and that derived by solving the classical 



vra 
aB 



equationlliP in the case of two defects. Here the distance between the impurities, kmp is 159a (a) 
and 79a (b) and l/{qQa) is chosen as 11.8 which corresponds to undoped CuGeOs. The solid curve 
expresses 9ci{x) = F{x) + tt/2 and the dotted curve is obtained from the classical equations with 
the boundary conditions of Oci{0) = Odikmp) = ^(0) +vr/2. Note that we have chosen the region of 
F{0) = F[limp) ~ —tt/2. Figure 2 shows the following facts. Our treatment that the cosine term 
is replaced by the term proportional to (j)^ underestimates the correlation between the impurities. 
Our result is valid for the case of the small doping, qoa/5 ^ 1. 
The free energy of eq.(^) per unit length, ftr = Ftr/L is easily calculated as 



2 PL ^ 

+ !^^2 _ B^^-\H(0,0;T)}_ f d^cOsF(x) 

a L J 

-^F(0,0;r), (15) 

where Eg = Va\J + is the excitation spectrum of spin fluctuation. In the following, we consider 
the case of T = and q^a/S » 1. Then the energy per unit length, etr = ftr{T = 0), can be 
calculated as follows, 

faqo\^ Vg aTT^ / aqo\ 2K ^ 



etr = - — o — + — o—d — + u 



27ra2 V 2 / tto^ 8a 12 



a 



1/2 
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Fig. 2. Spatial dependence of the classical part of the phase, 9ci, in the case of two impurities, where the distance 
between the two impurities are chosen as 159a(a) and 79a(b) with a being the lattice distance. The solid ( dotted 
) curves express the classical part of the phase obtained from the present study ( the classical equation in ref . 



where Cm(7r/2) = dt{l - cost)/t ~ 0.556811). In deriving eq.(16), we used qoa <^ 1. By 
minimizing eq.(16) with respect to u and ago/2, and defining and uq as the spin gap and 
the lattice dimerization in the absence of the disorder, quantities As/ and u/uq in the case of 
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J6/As <^ 1 are given as follows, 



A°5 




(17) 



Note that eq.(18) is consistent with eq.(14) in the sense that the spatial average of eq.(14) leads 
to eq.(18). Thus the spin gap and the lattice dimerization are suppressed linearly as a function of 
doping as is seen in the transition temperature. The ratio, As/Tgp, is given as 



Tsp Tl^y \8A^s C^Tl^J^j- (''^ 
For the case of Cui^s'Z'ia.sGeOs, the quantity As/Tsp decreases with increase of the doping rate, 
i.e., As/Tsp ^ A^/r5p(l — 175). The decrease of the spin gap originates from the fact that the 
coupling between the spins decreases due to the suppression of the lattice dimerization around the 
spin defects. Such a result indicates that the model having the constant lattice dimerization with 
the value in the absence of the impuritie^^ is hard to explain the suppression of the spin gap. 

§5. Summary and Discussion 

We have derived the action for the SP system in the presence of the spin defects, and have shown 
that the lattice dimerization is suppressed and the AF moments appear around the location of 
the defects. The doping dependences of both the transition temperature and the spin gap were 
calculated based on the action. 

The action was obtained from the Peierls-Hubbard model which includes not only spin but also 
charge degrees of freedom, because the substitution of non-magnetic ion for magnetic ion influences 
both degrees of freedom. The resulting action given in eq.(3) is equal to that of SP model with 
spin-lattice coupling which changes the sign at the location of the defects due to the kink of the 
charge. 

The transition temperature of the present calculation shows good agreement with the experi- 
mental result of Cui-^Zn^GeOs for 5 < 0.02. However there is a discrepancy in the sense that 
the saturation of Tsp appears in the experiments for 6 > 0.02 while Tsp of the present calculation 
decreases exponentially by the increase of 5. Further theoretical studies are needed in order to 
understand the saturation of the transition temperature of the SP state under doping. 

By applying the variational method to the action, the coexistence of the AF state and the SP one 
at r = was demonstrated without the assumption of the boundary condition at the spin defect. 
The present results are consistent with the assumption and the result by Fukuyama et al. . In 
addition, it was shown that the spin gap and the spatially averaged lattice dimerization decrease 
by the doping. The decrease of the spin gap is due to the suppression of the lattice dimerization 
near the spin defects. 
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Finally we comment on the low energy excitation in the spin gapQ'&ll2lll3EJ'E3). It is very compli- 
cated to discuss the excitation by the variational method. However the existence of the excitation 
in the present action of eq.(3) could be understood qualitatively by the following discussion. Ac- 
cording to the renormalization group analysis, the coefficient B in eq.(3) tends to strong coupling 
for K„ < 4, which shows that the system described by eq.(3) with Ka^ < 4 belongs to the same 
universality class. In the special case the Hamiltonian of the spin sector correspond- 

ing to eq.(3) can be mapped into that of non-interacting Dirac Fermions systems with the mass, 
mo(x), 

n = v^ J dx ^^ipl{-id^)Tpi - 7pl{-id3:)ip2] 

+ Jdx {mo(x) (^^1^2 + i^Ui) } , (20) 

where mQ{x) = 7raBm{x)u{x). The term proportional to cos2</)_(_ is neglected because the term is 
less divergent compared to that proportional to cos(/)+. In the case where mo{x) takes periodically 
two values (po and —(pQ with the same interval lint, the gap near the zero energy vanishes, because 

2/ 

mo{q = 0) = /q dxm{x) = 0. We also note another case that mo(x) takes two values (po and (pi 
alternatively with the random distribution of the interval length /. When the distribution is given 
by /o(0 = 6'(/)noe~"o' and = 6{l)nie~^^\ i.e., l/riQ and 1/ni being the mean length for 
and 01 respectively, the integrated density of states and the localization length can be calculated 
exactlyill' . Fabrizio and Melin@ have noticed that the Hamiltonian of SP system of XY chain with 
doping of Zn is given by eq.(20) and showed the midgap states by applying the above method. For 
the quantitative discussion, one must take into account the fact that the spatial variation of the 
mass, mo(x), is not exactly the telegraph type but changes the sign with the characteristic length 
of the order of l/q^ due to u{x) oc cosF(2;), and that the magnitude of the mass depends on the 
doping rate. The further exploration is needed in order to understand the low energy excitation in 
the presence of spin defects. 
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